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ABSTRACT 


The interaction of a single Lamb vortex with a free surface is analyzed 
numerically through the use of a finite-difference technique. The individual 
effects of gravity, viscosity, and surface tension are investigated within the 
range of the applicability of the phenomenon and the code used. The vortex 
is allowed to build up to its full strength in a relatively small time and then 
the evolution of the free surface, streamlines, and other details of the flow are 
calculated. The results have shown that the smaller the proximity of the 
vortex to the free surface, the larger the scar produced on its down-wash side. 
The effect of the surface tension is to reduce the amplitude of the free surface 
elevation. The viscous effects appear to be relatively small even though the 
calculations are, out of necessity, confined to a limited range of the governing 


parameters, as in all finite difference calculations. 
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I. INTRODUCTION 


The study of the flow about a submerged foil, in general, and a foil in the 
vicinity of a free surface, in particular, has a number of applications in 
hydrodynamics. For example, the effect of the free surface on the 
performance characteristics of partially submerged bodies, the stability of flow 
about control surfaces, and the surface signature created by the wake of a foil 
have emerged as important problems during the past decade. In these 
applications, the effects of gravity, viscosity and surface tension are almost 
equally important. The gravitational forces control the shape and the 
amplitude of the surface disturbances, viscous forces control in a more subtle 
manner the curvature of the surface and the rate of decay of the flow 
structures, and the surface tension controls the magnitude of the gradients 
and strains at the surface, and, thereby, the amplitude of the disturbances, 
together with the gravitational forces. It is, therefore, of paramount 
importance that all numerical calculations attempting to predict the behavior 
of flow resulting from a submerged or partially submerged foil account for all 
three effects in a three dimensional flow and, if not yet possible, in a 
relatively more manageable two-dimensional flow situation. 

As far as the recent naval hydrodynamic applications are concerned, the 
emergence of the remote control and observation technology, combined with 
the emergence of surface signatures, precipitated by the motion of near- 
surface submerged bodies, gave rise to the explorations of the non-acoustic 


detection of submerged bodies. Even though the technology is rather new 


and the exploration of its potential requires the understanding of many 
satellite technologies, the preliminary investigations have shown its 
potential and forged many interesting studies in the past few years. 

The first such investigation was undertaken by Sarpkaya and Henderson 
(1984) who have shown that the interaction of an ascending vortex pair with 
a free surface gives rise to scars and striations. The scars are relatively deep 
and long depressions on the downwash side of each vortex and are comprised 
of many vortices with axes normal to the free surface. These scars are bridged 
by the so-called striations which give the appearance of the rungs of a ladder. 
The character of these scars and striations are such that they give rise to 
coherent structures whose scale and frequency band are in the range of those 
of the Brag frequency. Thus, the alteration of the scale of the near-surface 
coherent structures (e.g., vortices or whirls) by single or double trailing 
vortices serves as a Brag-frequency discriminator and provides the connection 
between the non-acoustic detection of submerged bodies and the study of 
surface signatures. 

The realization of the foregoing led to considerable work in recent years 
on the interaction of vortices, wakes, and other types of flow with deformable 
fluid interfaces, all characterizable by the more general phenomenon of the 
interaction of vorticity with rigid and stratified deformable boundaries. 
Experimental work has ranged from the study of the interaction of submerged 
jets with a free surface, to the study of the interfacial turbulence in grid-stirred 
tanks. Numerical investigations have ranged from the use of vortex 


dynamics for inviscid flows to the use of finite-difference techniques for the 


solution of two-dimensional viscous and incompressible flows (with or 
without linearization at the free surface). 

The interaction of an ascending vortex pair with a free surface was 
analyzed in detail by Ohring and Lugt (1991) through the use of the full 
Navier-Stokes equations for an incompressible two-dimensional viscous 
flow, including surface tension and surface-generated vorticity. Out of 
necessity, their calculations were confined to relatively small Reynolds 
numbers (suitably defined in terms of the velocity of the ascending vortices , 
the initial vortex spacing and the fluid viscosity) in order to maintain stability 
and convergence. The effects of the Froude number (representing gravity) 
and those of the Weber number (representing surface tension) were properly 
accounted for. It is of importance to note that the value of the Reynolds 
number is not of particular concern for their investigation since their 
objectives were not so much the simulation of the ocean-truth or of the 
establishment of a practical predictive model, but rather the understanding of 
the physics of the phenomenon as to the size and evolution of the scars and 
the effect of the secondary vorticity on the motion of the primary vortices. 
The inclusion of the effects of surface tension were to shed further light on 
the anticipated effects of surface contamination which invariably exists in the 
oceans. 


Ohring and Lugt have shown that: 


High and low Froude numbers represent the two extremes of free surface 
yielding and stiffness, respectively. For an intermediate Froude number, 
a special rebounding due to the presence of secondary vortices has been 
observed: the path of the primary vortex centre portrays a complete loop 
(Ohring and Lugt, 1991, p.47). 


They have also shown that for relatively high Froude numbers (i.e., larger 
deformations) and for Reynolds numbers larger than about 50, the predictions 
do not appear to be strongly dependent on the Reynolds number, within the 
range of Reynolds numbers considered by them. In other words, the free- 
surface deformation without the influence of the surface tension is 
determined strongly by the gravitational effects or the Froude number. 
When the surface tension effects are included, the larger the Weber number, 
the larger the attenuation of the surface elevation, as shown by Ohring and 
Lugt (1991). This, however, raises the important question that the delineation 
of the effect of the Reynolds number on the results might depend not only on 
the magnitude of the Froude number but also on the Weber number since 
their increasing values produce opposing effects of the surface elevation. 
Thus, if the large Froude numbers give rise to larger elevation and thereby 
lessen the effect of the Reynolds number and the larger Weber numbers 
reduce the surface elevation and tend to intensify the importance of Reynolds 
number, it would be rather difficult to ascertain their respective roles when 
both the viscous and surface tension effects are included. Nevertheless, the 
choices for the numerical analysist dealing with the Navier-Stokes equations 
are limited partly by the stability and convergence of the solutions and partly, 
and perhaps more importantly, by the stability of the phenomenon itself. 

As far as the effect of a single vortex on the free surface is concerned, 
only Tyvand (1991) considered a vortex placed suddenly near a free surface 
and attempted to calculate the resulting wave phenomena for an inviscid 
case. His calculations, valid only for very short times, have shown that the 


sudden introduction of the vortex does not give rise to waves. This fact has 


already been established by others. Tyvand (1991) did not show any longtime 
calculations, streamlines, or surface deformations. 

The present investigation is a detailed solution of the two-dimensional 
Navier-Stokes equations for a viscous, incompressible fluid with a viscous 
vortex placed below the free surface. The finite-difference code used in the 


calculations will be described in the following sections. 


Il. NUMERICAL REPRESENTATION 


A. INTRODUCTION TO PROGRAM 

The numerical analysis discussed herein was performed using the 
RIPPLE computer program which was specifically designed to model 
"transient, two-dimensional, incompressible fluid flows with surface tension 
on free surfaces of general topology" (Kothe, Mjolsness, and Torrey, 1991). 
The code was modified to the extent necessary to solve the vortex/free-surface 
interaction problem and to produce the numerical and graphical results. All 
of the subroutines required for performing the numerical solution were 
organic to the program. 

RIPPLE solves either the Euler or the Navier-Stokes equations for 
incompressible fluids (with or without surface tension), using a two-step 
projection finite difference scheme which is second-order accurate in space 
and first-order accurate in time. It is capable of solving in either a Cartesian 
or cylindrical two-dimensional coordinate system, using a uniform or varied 
mesh, and, has the ability to model curved boundaries and internal obstacles 


using various types of boundary conditions. 


B. DESCRIPTION OF PROGRAM 

The following brief description of the numerical method follows closely 
that given by Kothe, Mjolsness, and Torrey (1991). 

The governing equations are the continuity equation for an 


incompressible fluid, 


(1) 


and the momentum equations, 


aV —- i 1 eae = 
—+Ve(VV)=-—Vp+-—Vet+g+— (2) 
= niell.o (NN) mR. B+ oF 


where 7, is the viscous stress tensor for a Newtonian fluid, p the fluid density, 
p the pressure, g the acceleration due to gravity and Fp the total body force. 
The particular properties of the fluid are assigned positions within a given 


computational cell as shown in Figure 1. 


Vi j+1/2 


Zz 
j+1/2 


Vi+/2,j 





y i+1/2 


Figure. 1 Cell Definition in the Grid (Kothe, Mjolsness, and Torrey, 1991) 


For the two-step projection method, Equation (2) is discritized with 


respect to time as, 


n 1.,ntl, 1 ne Sie 
— “on += vaele woe toned (3) 


The only implicit term in this equation is pressure—all others terms are 
approximated according to their values at the previous time step. The first 
step of the projection scheme solves for a velocity field based on incremental 
changes in the explicit terms of Equation (3). The second combines this 
intermediate velocity field and the pressure term of Equation (3) with 
Equation (1) to calculate the velocity field at the next time step. This 
combination forms a Poisson’s equation which is solved using an incomplete 
Cholesky conjugate gradient (ICCG) matrix solver. The momentum 
advection term of Equation (3) is solved using the second-order upwind 
method of van Leer. 

When solving for the free surface, the program assumes that the surface 
tension, 6, is constant and that the viscous stresses at the free surface may be 
neglected. This allows for the application of the continuum surface force 
(CSF) which represents the surface tension as an equivalent volume force 
rather than a pressure jump (Brackbill, Kothe, and Zemach, 1992). Thus, the 
free surface discontinuity is replaced with a smooth transition, or “color,” 
which varies across any computational mesh cell containing a free surface. 
The volume force is non zero only within these cells and is solved as one 
component of the body force term in Equation (3). A contact angle, 8, may be 
defined at any boundary and used to approximate wall adhesion forces. 
Because surface tension is solved explicitly, it is subjected to a linear stability 


time step constraint. 


Numerous techniques have been developed to perform the non-trivial 
task of modeling the free surface of a fluid. The method employed in RIPPLE 
is the Volume-of-Fluid (VOF) method. VOF is a donor-acceptor differencing 
technique which makes use of a scalar field function, F, to act as a marker of 
the free surface. F describes the fraction of a cell which is occupied by the 
fluid: 

F(x,t) = 1 in the fluid; 
F(x,t) = 0 in the void; and, 


0 <F(x,t)<1 at the free surface. 


The free surface, therefore, is not continuous, but rather, is represented by a 
series of discontinuous line segments described by the value F in the surface 
cells. Although, a reconstructed free surface is not required for surface 
tension calculations (since this is handled by the CSF model), it is needed to 
ensure that the fluid in the vicinity of the free surface is accurately advected 
over the domain. Both the reconstruction and advection of the VOF function 
is performed using a Hirt-Nichols (H-N) algorithm. Because the H-N 
algorithm is explicit in time, the maximum time step is subject to the 
Courant condition. Also, the fluid color, used in the CSF model is related to 
the VOF function. The viscous stress term is solved explicitly in time using a 
backward difference scheme, so it too is subjected to a linear stability time step 


constraint. 


C PROGRAM PARAMETERS 

The first parameter which must be chosen is the size of the domain. 
Because it is desirable to limit the problem run-time to some reasonable 
length, excessively large domains are impractical. Therefore, the size of the 
domain must be limited, while simultaneously minimizing the impact of the 
resulting boundaries which must be in relatively close proximity to the 
region of concern. Given an initial vortex depth below the free surface of ho, 
it was found to be sufficient to place the side boundaries at t4hy and the lower 
boundary at 3h, below the vortex center, for a Froude number Fr (= 
T/hg4/gho ) = 7.5. For Fr = 13.8, the side boundaries were placed at +6hp and 
the lower boundary at 5hy9. The top boundary needed to be located high 
enough above the free surface so as not to interfere with the largest vortex- 
induced free-surface deformation. For this analysis, 1.7hg was deemed 
sufficient. 

Closely related to the size of the domain is the type of boundary 
conditions imposed on its periphery. RIPPLE is capable of applying free-slip, 
no-slip, periodic, applied pressure, continuative outflow or a specified 
inflow/outflow condition to any of the four boundaries. The last three 
conditions have no relation to this problem. The no-slip condition is not 
appropriate because it enhances the influence the boundaries have upon the 
domain. The periodic condition simply reflects the velocities, pressures and 
free surface positions about the boundary and does not accurately represent 
the images produced by a single vortex. The free-slip condition is ultimately 
the best choice since it limits the effect of the boundaries and crudely models 


the vortex images on the sides and bottom of the domain. 


10 


The grid size can be either uniform or varied in either direction. 
Because a varied grid is applied along the entire length of the domain (for 
example, a variation in the grid in the y-direction must be applied from the 
bottom to the top of the domain) the number of cells are increased along the 
entire strip and not just in the region desired. It is as a consequence of the 
decrease of the cell size in such regions that the demand for CPU time 
increases. Thus, a choice must be made between the increase of the CPU time 
and the increase of the accuracy of the computation over a “cross-shaped" 
refined region. A uniform mesh was, therefore, used with its coarseness 
limited by RIPPLE’s ability to accurately solve for the free surface for the given 
values of the Froude number (Fr = T/hoVgho ), the Weber number (We = 
o/ pgh?), and the Reynolds number (Re = hovgho /v). A 192x112 mesh grid 
was capable of meeting the needs of this analysis within a reasonable run- 
time using a SUN SPARCstation 10™ (25 to 70 hours). 

The program provides for a choice between a fixed time step and a 
variable time step (adjusted automatically during the course of the 
calculation). Regardless of which method is chosen, the code requires that an 
initial time step and a maximum allowed time step be provided. For this 
analysis, the variable time step was used because it reduced the run-time 
while making sure that the linear stability requirements were met. It was also 
determined that an initial time step approximately equal to 1/10-th of the 
maximum time step of the three time constraints was sufficient. For initial 
time steps larger than that, the code frequently exceeded one of the time 


constraints, requiring it to recalculate the current time step. 
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As was the case for the grid size, the tradeoff for defining the required 
convergence criteria for the Poisson pressure solver was the solution-accuracy 
versus problem run-time. Through limited trial-and-error, it was found that 
a value of 0.001 was acceptable. 

As described previously, the advection calculation is performed explicitly 
with respect to time and must, therefore, meet the Courant criteria. Implicit 
in the program is the assumption that all time steps have Courant numbers 
less than 0.5 — any step not meeting this criteria must be recalculated. To 
prevent an inordinate amount of recalculation, the Courant number was set 
to 0.38. This ensured that the Courant number after the calculation was less 
than 0.5 and that the advection solution remained stable. The program also 
allows for choosing the type of momentum advection. Based on the 


recommendation of the authors, van Leer’s method was used. 


D. ADDITIONAL PARAMETERS 

A code such as RIPPLE allows the user to define numerous parameters 
and variables, but it obviously cannot anticipate every particular situation to 
be encountered. The interaction of the vortex with the free surface had to be 
correctly implemented and incorporated into the code, together with the 
selection of the appropriate boundary conditions and parameters which are 
compatible with the demands of the flow field. Therefore, it was necessary to 
include a routine which defined the desired vortex and the variables 
necessary to vary the governing parameters associated with the simulation 


(Fr, We, Re and @=ro/Nhg, the ratio of the core size to the depth of the vortex). 
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A subroutine designed to specify the boundary conditions was modified 
to characterize a single vortex. The swirling motion was prescribed over a 
square region, surrounding the vortex core. The velocities were solved using 
the velocity equation derived from the Lamb vortex model in Cartesian 
coordinates, 


esi 
Mt) | 1 expd—1.2544] Y = (4) 


geht) = — 
2ny y* +27 Ig 


where I is the vortex strength, ro the vortex core radius and y and z as shown 





in Figure (2). 


Vy,z,t) 


MA 


Figure. 2 Definition of Vortex Characteristics 


The v and w components of velocity are then given by, 


v(y,z,t) = V(y,z,t)} ——_—— (5) 
y2 + Z 
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and 


w(y,Z,t) = =Wy 279 eer (6) 
y + Z 


This required the ability to define I, ro and the location of the vortex center, 
(Yc, Zc). The length of the sides of the square region in which these velocities 
were defined was set equal to 3rpo, with 1/3 <a<1/2. 

It was recognized that to place the vortex instantly at full strength at time 
T (= ho/ Ve ) = 0 would likely "shock" the system, thereby requiring long run- 
times to allow an equilibrium condition to establish at the free surface. Thus, 
the vortex was initialized at some fraction of its ultimate strength, and 
allowed to build up linearly over a predetermined time period. Typically, the 
vortex was initialized at 0.01I° and allowed 7T to 10T normalized times, to 
reach full strength. It was desired that the vortex/free-surface interaction be 
examined in a "natural state," where not only was the free surface shaped by 
the vortex, but the vortex distorted and moved by the free surface. It was, 
therefore, necessary to cease providing the vortex with a source of circulation 
so that both the vortex and the surrounding flow field could reach a non- 
forced state of relaxation. Subsequently, the quasi-steady state of the 
interaction was examined. To achieve this, one last parameter was added 
which specified the time at which the velocities within the square region 
surrounding the core would no longer be redefined. For all of this analysis, 
this time occurred at 25T, as this was sufficient to allow the free surface to 


reach a state of dynamic equilibrium, yet early enough to prevent the 
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boundaries from influencing the interaction. After reviewing the graphical 
output from numerous runs, it was determined that the best time at which to 
examine and compare the above interaction was 27T. Most of the runs were 
continued far beyond the time 27T. An example of the evolution and demise 


of a single vortex at larger times is presented in the discussion of results. 
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I. PRESENTATION AND DISCUSSION OF RESULTS 


A. INTRODUCTION 

Among the numerous calculations performed, only six representative 
single-vortex/free-surface interactions will be discussed in some detail. The 
six runs are classified according to the combination of the Froude, Reynolds, 


and Weber numbers: 


TABLE 1: THE RANGE OF THE GOVERNING PARAMETERS 


Set #1 Fr = 7.50 Re = 40 We = 0.033 
Set #2 Fr = 7.50 Re = 20 We = 0.033 
Set #3 Fr = 7.50 Re = 40 We = 0.333 
Set #4 Fr = 7.50 Re = 20 We = 0.333 
Set #5 Fr=13.80 Re=40 We = 0.333 
Set #6 Fr=13.80 Re= 20 We = 0.333 


Each set of calculations show the streamlines, velocity vector plots, 
v/ Vmax versus y/No, W/ Vmax versus y/No, V/ VMax Versus Z/No, and w/VmMax 


versus Z/No. 


B. ASSESSMENT OF THE EFFECTS OF VISCOSITY 
As noted from the foregoing table, the first two sets of runs are designed 
to evaluate the effect of the Reynolds number, however small the range of 


the Reynolds number may be. The comparison of Figures 3 and 9, 4 and 10, 5 
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and 11, and 6 and 12 show that aside from the secondary details in the vicinity 
of the vortex core, the shape and magnitude of the free-surface deformation 
are nearly identical. The secondary differences, attributable to the rate of 
diffusion, show in the straining of the vortex, in the vorticity contours 
around the primary vortex, and in the shape of the streamlines surrounding 
the primary core. It is noted that for the case of the smaller Reynolds number 
(Figure 9), the lower part of the main vortex is comprised of two cells of 
streamlines whereas the one for Re = 40 (Figure 3) is less diffused and the 
streamlines below extend over larger regions. In other words, there are subtle 
differences between the two streamlines in spite of the relatively small 
difference between the two Reynolds numbers. It will be shown later that the 
said differences are further enhanced with increasing Froude number and 
decreasing Weber number. 

Figures 5 and 11 show the v-component of the normalized velocity on a 
transverse line passing half-way between the original positions of the free 
surface and the vortex center (z/hNp = 0.50). The magnitude and direction of 
this velocity along the y-axis depends on the particular position of the z = 
constant plane. For a non-deformaing free surface, it can be demonstrated 
through the use of a Kelvin oval that the velocity should be maximum 
directly above the vortex and decrease symmetrically as the distance from the 
vortex increases as shown in the insets of Figures 5 and 11. The reasons for 
the change of direction as well as the difference in the magnitudes of the 
negative values of the v-component are a consequence of the asymmetric 
deformation of the free surface on either side of the z-axis. The rise in water 


elevation and the drop in the v-velocity to negative values on the upwash 
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side of the vortex are both larger for reasons which can be easily explained by 
examining the streamlines and the velocity vector plots. 

The w/Vmax versus y/No is shown in Figures 6 and 12 for the governing 
parameters of the first two sets of calculations, shown in Table 1. In these 
plots several facts are evident: (i) the velocity profile is not odd-symmetric 
with respect to the z-axis, as it would have been had the free surface remained 
undeformed; (ii) the magnitudes of the regions to the left and the right of the 
z-axis are significantly different, again due to the asymmetric deformation of 
the free surface; (ii) for the smaller Reynolds number case (Figure 12), the w 
profile is relatively more symmetric, as would be expected; and (iv) in both 
cases, wW approaches zero from negative values on the left and from positive 
values on the right. 

The transverse component of the velocity, v/VmMax, versus z/NRg is 
shown in Figures 7 and 13 for the two calculations under consideration. 
Aside from their magnitudes, the two profiles are virtually identical. It is 
noted that the vortex center has shifted slightly downward in time due to the 
mutual interaction. It is also noted that the velocity profile is not symmetric 
with respect to the z = 0 line because of the proximity of the free-surface. 
Furthermore, the ratio of the velocity extrema is larger than unit for the 
larger Reynolds number case whereas it is smaller than unity for the lower Re 
case. 

Finally, the variation of W/VmMax with z/ho is shown in Figures 8 and 14. 
Normally, this velocity should be zero for a non-deforming surface. Thus, 
however small, the magnitude of w/Vmax is one of the most important 


measures of the deformation of the free surface, i.e., the larger and more 
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asymmetric the deformation is, the larger the deviation of the w-component 
of velocity from zero. Figures 8 and 14 show that for a vortex rotating 
clockwise below the free surface, w becomes plus below the vortex axis and 
negative above the axis. In other words, there is mass flux into the vortex 
core in the z-direction. It can be shown, however, that the net mass flux into 
the core around a circumference is zero by virtue of the solenoidality of the 


flow. 


C ASSESSMENT OF THE EFFECTS OF SURFACE TENSION 

Now, the calculation sets of 3 and 4 (see Table 1) will be discussed. It 
should be noted that these two sets differ from the first two only in terms of 
the Weber number. Thus, the sets 1 and 3 and sets 2 and 4 have identical 
Froude and Reynolds number and any differences in their characteristics will 
be indicative of the increase in the surface tension or of the Weber number. 
In assessing the role of the Weber number on the flow behavior, it should 
also be noted that the surface tension has been increased by ten fold whereas 
the Reynolds number was only doubled. 

Figures 15 and 21 show the streamlines for cases 3 and 4. The first 
striking difference between the two figures is the magnitude of the 
deformation of the free surface: the smaller the Weber number, the larger is 
the deformation. This is somewhat expected on the grounds that the 
additional surface tension gives rise to surface tension gradients and surface 
vorticity and these in turn absorb energy which would have been otherwise 
stored in the fluid as potential energy in raising the free surface. A similar 


conclusion can be drawn regarding the low Reynolds number case, as seen in 


iy 


Figures 9 and 21. A closer examination of the streamlines very near the 
region where the free-surface slope changes (nearly above the vortex) show 
that for small Reynolds numbers (Re = 20), the build-up of surface vorticity is 
somewhat stronger than the higher Reynolds number cases. This is expected 
on the grounds that the flatness of the surface and the diffusion of vorticity in 
the lateral directions should lead to secondary circulations between the vortex 
and the free surface. 

The vector plots of the velocity are shown in Figures 4 (Fr = 7.50, Re = 40, 
We = 0.033) and 16 (Fr = 7.50, Re = 40, We = 0.333) and in Figures 10 (Fr = 7.50, 
Re = 20, We = 0.033) and 22 (Fr = 7.50, Re = 20, We = 0.333). These figures 
show that even though such vector plots are valuable in pointing out the 
direction and magnitude of the velocities, they do not allow one to develop 
an integrated mental image to draw conclusions regarding the interaction 
between the vortex and the free surface. These can only be achieved through 
a careful perusal of the velocity distributions along selected lines and 
directions. 

The second most significant impact of the surface tension must surely be 
associated with and be a consequence of the flattening of the free surface. 
Namely, the velocities in the transverse direction are enhanced as the 
streamlines are confined to narrower regions above the vortex (e.g., at z/hp = 
0.5, as in the present case) and the v/ Vmax values should be larger for the case 
of We = 0.333 than for We = 0.033, as seen in Figures 17 and 5, respectively. 
Figures 11 and 23 exhibit the same behavior, but less dramatically. 

The vertical component of velocity along a transverse plane (z/No = 0.50) 


for the case of Fr = 7.50, Re = 40, and We = 0.333 (Figure 18) may now be 
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compared with that for Fr = 7.50, Re = 40, and We = 0.033 (Figure 6). 
Apparently, due to the elevation of the transverse plane the difference in w is 
not as accentuated as the v-component. In fact, the solenoidality of the flow 
for a given vortex position necessitates that the increased v-velocities (the 
bunching of the streamlines) should be compensated for by a smaller 
difference in the w component. For the case of the lower Reynolds number, 
Figures 12 (Fr = 7.50, Re = 20, and We = 0.033) and 24 (Fr = 7.50, Re = 20, and 
We = 0.333) support the same conclusion. 

A comparison of Figures 18 through 20 for (Fr = 7.50, Re = 40, and We = 
0.333) with Figures 24 through 26 for (Fr = 7.50, Re = 20, and We = 0.333) 
shows that the dependence on the Reynolds number of the particular 
velocities are not significant and certainly no more than that described earlier. 
It suffices to note that the most important difference is in the v-component of 
the velocity and the degree of the differences on all velocities, including v, 
depend strongly on the proximity of the transverse plane to the deformed free 
surface. Here, it is implied that the difference is only in the magnitude of the 
velocities and not in the physics of the phenomenon. In other words, no new 
physical events are expected to occur as one approaches the free surface. It is 
also important to be reminded of the fact that these calculations deal only 
with laminar flows. In turbulent flows, the interaction of vorticity with the 
free surface gives rise to coherent flow structures which cannot be predicted 
here through laminar flow solutions (Sarpkaya and Suthon, 1991). 
Nevertheless, the overall behavior of the free surface, the generation of free- 
surface vorticity, mean-flow characteristics, the regions of formation of the 


primary and secondary vortices, the relative influence of various parameters 
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and several other experiment-guiding results can be predicted through the 
use of laminar flow calculations. It is on the basis of such results that 
turbulent flow measurements can be made to delineate the character of 
coherent structures and the transfer of vorticity from the vortex to the free 
surface and the generation of new structures and their redistribution in the 


transverse plane in the form of quasi-coherent quanta. 


D. ASSESSMENT OF THE EFFECTS GRAVITY 

The next issue to be taken up is the effect of the Froude number on the 
characteristics of the flow. The calculations are presented in sets 5 and 6 (see 
Table 1). The role played by the Froude number (signifying the effects of 
gravity) will be illustrated through the comparison of the results between the 
sets 3 (Fr = 7.50, Re = 40, We = 0.333) and 5 (Fr = 13.80, Re = 40, We = 0.333), 
and sets 4 (Fr = 7.50, Re = 20, We = 0.333) and 6 (Fr = 13.80, Re = 20, We = 
0.333). 

Figures 27 and 15 show that there are two fundamental differences 
between them. First, the surface elevation is larger for the larger Froude 
number. Second, the secondary vorticity is much more pronounced. Both of 
these effects are somewhat intuitive and could have been predicted in gross 
terms. However, the exact shape of the free surface, the deformation of the 
vortex and the relatively larger excursion of the vortex center seen with the 
larger Froude number could not have been anticipated. Figures 27 and 28 
show that the vortex center, which was at (0, 0) at time zero, has gradually 
shifted to the third quadrant of the coordinate axis, but near the origin. 


Additional comments may be made regarding the streamlines and the 


vorticity distribution in the regions below the vortex. Even though they 
further show the effects of the magnitude of the Froude number, here the 
attention will be concentrated on the near surface structures. 

Figures 29 and 17 show that differences between the two Froude 
numbers is accentuated. Likewise, Figures 30 and 18 show that the overall 
character of the variation of the w-component is quite similar for the two 
Froude number cases, except for the fact the growth of the larger secondary 
vorticity, particularly to the right of the scar, gives rise to a ‘hump’ in the w- 
velocity as shown in Figure 30. This is expected to be further amplified in the 
transverse planes closer to the free surface. Figures 31 and 19 show that the v- 
component of the velocity is dramatically effected by the magnitude of 
deformation of the free surface as exemplified by the strong asymmetry of the 
top and bottom halves of the profile and by the absolute values of its 
magnitude. Finally, the w-component of velocity along the z axis is 
compared in Figures 32 and 20. As noted earlier, this velocity component is 
rather small for obvious reasons. Thus, its distribution is strongly influenced 
by any parameter, large or small. In this case, the migration of the center of 
the vortex, induced by the proximity and the degree of deformation of the free 
surface, could and has changed the distribution of the w-component. The 
vortex center in Figure 32 is such that this velocity component is negative 
through the entire z-axis, unlike the previous ones. The unidirectionality of 
the w-component along a single z-axis does not change the fact the mass flux 
through the vortex is maintained. 

It is clear from the foregoing that the Reynolds number effects are subtle 


and the degree of diffusion of vorticity exhibits itself in all velocity profiles. 
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The effect of the Froude number is to amplify the magnitude of the 
deformations, enhance the secondary vorticity generation indirectly through 
the increase of the radius of curvature of the free surface, and through the 
straining of the vortex. The effect of the Weber number, however, is opposite 
to its magnitude in practically all aspects of the evolution of the flow. The 
larger the Weber number, the smaller are the deformations. In that sense, the 
surface tension and gravity effects act in opposite directions. Thus, one would 
expect that the cases of relatively low Froude and Weber numbers (e.g., 
Figures 3-8) and the cases of relatively high Froude number and Weber 
number (e.g., Figures 27-32) should be comparable. In fact Figures 3 and 27 
show that the free surface deformations are comparable. They also show that 
the near doubling of the Froude number has a greater influence on the 
magnitude of the free surface deformation then the ten fold increase of the 


Weber number. 


E. GROWTH AND DECAY OF THE VORTEX 

In the foregoing, the effects of viscosity, surface tension, and gravity on 
the interaction of a vortex with a deformable free surface have been discussed 
at a specific time without regard to the transient state. These have been 
instructive in clarifying the combined as well as separate effects of the 
fundamental physical parameters and in paving the path to the design of 
physical experiments. Clearly, it would be equally important to know as to 
how a particular state has been reached and how the vortex decays with time 
if its circulation were no longer maintained constant. This section deals with 


these issues. 
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In view of the fact that it will be nearly impossible to discuss the 
evolution of a fluid state for all possible combinations of the governing 
parameters, only the most representative values of the Froude, Reynolds, and 
Weber numbers were chosen, i.e., Fr = 7.50, Re = 40, We = 0.033 (case 1). 

Figures 39 through 47 show, starting at T = 10, the evolution of the 
streamlines. At time T = 10, the vortex has already reached its full strength. 
The circulation is held constant between the times T = 7 to T = 25. For T > 25, 
the vorticity of the vortex is allowed to diffuse without the infusion of new 
vorticity into the core. Thus, the circulation in the original core is allowed to 
decrease with time. 

A cursory examination of the figures 39 through 43 (at T = 26) shows that 
the free surface quickly reaches its steady state shortly after T = 10. However, 
after T = 26 (Figure 43), the diffusion of the unreplenished vortex begins to 
manifest its effects on the entire flow in the form of decreasing surface 
elevation, uniformalizing of the streamlines (showing that the vorticity 
gradients are becoming less steep), and the scar becoming more smooth. The 
Surprising aspect of the decay process is that a steady state which has been 
essentially arrived at by T = 10 is largely maintained even at times as large as 
T = 30, at least in the immediate vicinity of the vortex. This is an indication 
of the longevity of the larger structures near the surface. The measurements 
previously carried out at Naval Postgraduate School (see, e.g., Neubert 1992) 
show that energy can cascade towards larger structures and thereby decay at a 
slower rate. That decay is certainly a function of a number of factors. Aside 
from the domain in which the vortex is born and subjected to decay, the 


condition of the free surface at a particular instant, the viscous, gravitational, 


AS, 


and the surface tension effects (surface contamination) will certainly 
influence the duration of the life of the vortex. If one were to assess the half- 
life of the vortex in terms of its decay to half the vortex strength, one might 
conjecture that this would be a very large time, not assessed in the present 


investigation due to CPU limitations. 
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IV. CONCLUSIONS 


The interaction of a single Lamb vortex with a free surface has been 
modeled using a finite-difference scheme. The values of the Reynolds, 
Weber and Froude numbers were varied and the results compared in an 
effort to assertain the influence of viscosity, surface tension and gravity upon 
this interaction. 

Compared to that of surface tension and gravity, the influence of 
viscosity on the evolution and steady-state of the surface deformation is 
rather secondary. Its major contribution is that of controlling the rate of 
diffusion and degree of vortex strain. 

Within the range of Reynolds numbers examined, the effect of surface 
tension is to moderate the effect of the remaining parameters upon the 
interaction. In most cases, the result of an increase in the value of the surface 
tension, or Weber number, is opposite to that of the Froude number. 

The gravitational force has the most notable influence on the 
vortex/free-surface interaction since it changes the curvature of the free 
surface with increasing circulation or with decreasing proximity to the free 
surface. Furthermore, the gravitational force indirectly controls the amount 
of vorticity generated at the free surface. 

All of the foregoing is influenced to varying degrees of intensity, as well 
as complexity, by the non-symmetric nature of the free surface. This leads, 
among other things, to the stretching and tilting of the vortex, the flow of 


surface vorticity to the downwash side of the vortex, and, most importantly, 


ay 


to the distributions of the v and w-components of the velocity. The v- 
component increases with the decrease of the vortex depth and with the 
flattening of the free surface. The w-component may be both positive or 
positive-negative depending on the elevation of the vertical line along which 
the velocity is calculated, the time of the calculation, and on the magnitude of 
the governing parameters. 

The results have demonstrated that the deformation of the free surface 
for two-dimensional laminar states can be calculated through the use of a 
finite-difference code, capable of including all three effects: viscosity, surface 
tension and gravity. 

It remains to be seen as to what the effects of calculation domain size, the 
character of the initially imbedded vortex (e.g., Lamb vortex versus a Rankine 
or Rosenhead vortex), the individual or combined effects of integration time, 
grid-size, and the boundary conditions be, among the most important 


computational input parameters. 


28 


APPENDIX 


(27 = L1® €€0'0 


IM Pure Op 


Or 


‘0S'Z 


= IJ) SaUT[WIRaI}S TH Jas 


Ou/A 


‘€ AINSI 





29 


(ZZ = L1e €€0'0 = 9M pur OF = ay ‘OG'Z = 4) S10}DaA AWOOJAA 1# 13S 


OU/A 


aaa 
a4 
a 


aase ue bB_u_U_aa 
aaa 
aa 
aa 


eye » ds 6 dB db DS 
ea db &6 & 3d 6 Sh 
>» am» &@ &6 & 6 6b 8 
aéaés66 8s 
aaéssbbi’b 
aasseb bib 
aéeé686464686 46464648 
aaaéaa ld 
a¢4é4 4644 46 46 
4¢444 4 8 6 
¢ ov @ ¢ € 4 4 4 «@ 6 


a@a@es24 939 4€93 4 «@ 


a @# a2nt 2 2 @ ea a2 «© ew wy 
a> ea ea #4 @e wt eaasevwrvwsewsewessia 
aa 


eaten ee 


,yrPPpererePrP Per er ere eyeereessrssssb db 


Pomme ww 4 4 4 fF 


DP iW. 


i A a a i 2 ee a | 
| A A a ee a er ek ee VY 
| i a A a a a ee ee We | 
J a ee a a ee es ee er Ve 
,rpPrPPrF er em mC PrhUh Oe Atroeewawvywvyerspasb’ 
i ee ee a er es er ee We | 

? 

qe0ea48 4&4 & @weewedevdev ti ¢€ 46 bh bebe 
q¢ea844 &4 @> PPwHwrevwedeettie)iteth*_—eoU_egas 


ft 4 me ewe erwreneaeee 4 & F 


- 

Ld 

¢v 

¢ 
¢ 
¢ 
¢ 
gv 

- 
bd 
= 


> P&H & He Oe 6 4H OH fl 


~ 
& 
~ 
ad 
& 
bal 
- 
& 
Lad 
~ 
Ld 


rrrevrt ep eo > ® 


"TeeawrtvereF Ff 


q 
t 
wv aga eee Ra wewevwewvevewveewvt ee bhbeeoe 


7 ¥44 4 eee mee eweeweverveveve et tee 


*#e#qroaett= 8 & &@ wpe vet ¢€ @¢ 4 eb se 


1°. P & eo Hh oo OO 


t 
t 


' @(@ Pe ae H8a ee oD ww & @& fb lll 
g 
t 


a 
ry 


7 aea 4 me eae mewmereereerevsevwsecdetetiéa 


vet 444 ema meereereervrerwververwvv et ets 


‘Pp amn3iy 


4 *@* & ea R&R SB wPererererwreerereeeveysrg¢ 
=> RH @&®» ew PwP Pe FP PP PP Pw 
a ewe PB PB PP FP PP FP FP PF PP PP Pe 


7a et 4 & me emenrreereveevewsevrdet et_"! 
T¢4¢e@eaeearnrperwreereerevvwvr¢«vs¢ 





30 


(ZZ = LJ" €€0'0 = 9M Pure Op = ay ‘OG'Z = 14) 


0G°0 = Sy /Z je adueysiq] asJaAsuesy snsiaA XV A /A [# JaS 
Ou/A 
OV O'? 00 OC Oi — 


‘¢ ain3Biy 





31 


(ZZ = Lie €€0'0 = 29M pur Op = ay ‘0G'Z = 1) 


0S°0 = Oy/Z ye adue}sIq aSJaASUBIT sNSIaA XVNA /M [# 19S 


u/A 
O’v O°? 00 Oe O'v— 


‘9 AINSI] 





32 


(ZZ = L1" €£€0'0 = 9M Pure OF = oy “0G°Z = AH) 
00'0 = CU/A Je aduUR\sIG JedNIVA SNSIAAXVYNA /A THIS =*-Z aINBIY 





(ZZ = .L 38 €€0'0 = 9M pur Op = ay ‘0G'Z = 1) 


00°0 = °Y/A je aouRysSIG JeoVaA SNsIaAXVYNA/MTH#YAS = -g aNBIy 


XVNA/M 
L°O 0'O L°O- 





(ZZ = L318 €€0°0 = 
L 1 €€0'0 = 9M pure 07 = dy ‘0G'Z = 4) SaUTTWIRAI]S ZH 9S «= ‘6 AINSI 


OU/A 





35 


(ZZ = Li €€0'0 = 89M pur 0Z = ay ‘0G'Z = 4) SIOJDIA AWOO[AA 7H IOS ‘OL ainBiy 


O'V 


e 
4 
4 
e 
e 
i 
e 
A] 
| 
e 
v 
e 
A] 
e 


a a2 a2 42am 24282 @2 428 2908928424 


a aaeaeanaaaasaf4£284 84 4 


eo #@#@ 4482 424824244 UY & 


QO # #@ #4 #4 HA 2A PA wm maoeanoyvwewvsseroiéia 


oe + bP # ae #@ A wee mae eaeaaswsewsyreedss sh 


| a a a ee ee ee rr Y 


| a a a ee ee a | 


| a a a a a ee es ae Wr | 


| J 2 a da 2 2 ee ee ee we | 


| At 2 2 2 2 a 2 ee” 2 es ee fe | 
PP PrP PrPeeremhlarmhlc hhh Ph PO Aa rFM 7a owseserhrdssaa isa 
PPP PrP rem rr oO he eo owsweeyrd='h bh db da ib 


PPP PrPer ttt P 
> Ppt? € 24 2 FE 


a2ererv»espéaéddbésdé 


es & & & Od O88 


»ep dBd dB &’ & 6 
aanraéseaeab db 


PF pPrFeee @eaeaese 4 


> tT & me Pee eeese & 4 


- rr Fe eea ea & & & 


CU/A 
0°O 


ao@wst 9 4a @ aw @ 
w-s8t wsgswaass990 es 
w-wt ws sas wa vw 4 @ 
bt in ee i | 


swsewts3s3t3tbshGhCU CU 


aos Fe ewe wa ws 


a4a4a644 4 


arv%* Ste eeer sense er * 


VaoaVewreeweaesp sh & & 


eer Wee ee a 
oeeeuanaa 


4éaée«¢686866 6 44 


tt &@ #4 Ha hh OG 4h eh 6 hf 4 fo 


ee ee ee ee ee ee ee ee ee | 


@¢¢€ @ PO PRhlPMrmhlClMelUlUlUD 
@e#0u¢? @¢ Pp a 
@ vegeta? € 4 


eo ¢€¢6 64644644 4 6 6 


@*¢eet¢ 44 & @ mw wp wpgegegeeet besebb_b_('t 
7 @qeee¢44 4&4 &R mewmpyeevg¢e¢eeteéethshebeaet 
T@et41wtete eae rwmernempmevwvvreveerertibaa 


¢€¢¢4484 08 82 PP ePwesvewewewe ht bb tb 


@#¢@€ 44646480882 Pw eww wv vvev 4£ ees b 


*¢ 44 eae &e meme wrPerreerveervewv eee 


*¢¢4464 808 eae PrP ew rw we ww wwe ee fb 


i 2 2 | 


Ov 


a &*® & @a@eaoerrerrerervewFevg¢ 


a &®# &# ®&# Pe RP ewrPrreerperePr we 


= PrP PrP Pe Pere PrP PP FPrhlUlcrlUr hh 





Do 4 


36 


(ZZ = L 38 €€0'0 = 9M Pure 0Z = 2Y ‘OG'Z = A) 


0S'0 = CU /Z }e adURIsSIG aSIBASURIT SNSIAA XVN A /A THYBS TT ans] 


Cu/K 
O'Vv O°? 0°O OrG= 0'v— 





37 


(ZZ = L 18 €€0'0 = 9M pue OZ = AY ‘0G'Z = 4H) 


0S'0 = °Y/Z Je aduR\sSIG asIaASURIT SNSIAA XVNA /M ZH 49S ‘ZI an31y 


Ou/A 
O'v O°? 0'O OA O'v- 





(ZZ = L 1" €€0'0 = 89M pur OZ = ay “OG'Z = A) 


00'0 = °U/A Je aduR}SIG [eITIAA SNSIDA XYN A /A ZHYIS “ET AINSI 


XVNA/A 
G0 0°0 G 0- 





39 


(ZZ = L3€ €£0'0 = 9M Pue 0Z = AY ‘0G°Z = 14) 


00'0 = Cy /A ye aouLysiq JeondaA snssaA XYNA /M ZH 9S PL aINBLY 





40 


(ZZ = L 38 €€E'0 = 9M PUR Of = ay ‘0G'Z = AY) SaUTTUTVANS C4 }eS “ST aINSIy 


°u/A 


7 / oF : ws vel 
yt e ‘, ® 
AH) ogy 
é eds Po --— ee 
Hh + LG S\\ i i‘ 
4! “ue rif Nee 3. 
Hae ee eu Py j .* 
bus 14 18 ii wb 
is vy Apr i 
int Nn Ss Ae 
‘ Ae > 3 e 
3 * s 5 
*. ' a 
AS ; « . 





4} 


(ZZ = L 1 €€E0 = 9M PUe OF = ay “OGL = AY) SIOPAA AyDOJaA E# 9S ‘OT AINBHT 


OU/A 


7~ yo debe bh db 8 
we» 6 6 &6 Ob DS 
aaaeaubat’ 
aseeaeebeaua’ 
aaebeeseb ab 
aebee@e asd 
a¢é44644686 6 6 4 
ad¢éé 4446 64 6 4 4 
4¢4664444686464668 
4é@aé 


-~e vvoriba b&b 8 
q 

g 

q 


q¢¢@e¢ets@#é@#é@éesea 


aaqaaeuwwe¢@s 
aa 4aae37+7Fwweesey 
7¢a@e# 44 4& ©&e@ wo wpevvwsedsed¢ ttt tba 


a 


a ewfeerFre PrP FP FP FP FP FP PF FP PP 


—~_> 


— 
Pad 


i a ee 2 | 
i a a ee a ee ke Wr a 
& ma m@e @eereeeevevrevg¢g 

a @&# @&# Fw weeereeFerereeree pes 


Ow Pomem~ea~w en a 4 4 tf 


4 ¢«¢ a & me e@eweewreeweeweseeseeev¢e¢ 


\ 
4 
¢ 
¢ 
¢ 
4 
4 
4 
a 
4 


A A a A a i ee ee Ve 
rer FF mena e__ @ ft 1 


>rPPrPPrPPerFPrPrP eH FA HA Het BDewvweyerdseab a 
A A A a a a a a 
yrPPrPPrpPererr PrP Pr heh hh re reece erserwsadba b 

PP Aanaaeeantwtetaas a 

qe<@4*4¢€4*® @eevwprevedwee¢et#t@#ées_she bbb 
1 @¢446444%4 @H@# weovrvwewewcte¢e¢ecettbtibe_iaie_ bp 


t* € 44 82 S&S BO wPwMPwrVwewrerewwweee ea’ 
© 444 &@& me eae ww we Kw wPwrrwrwwewee? 


a 
an 
a + 
r" ¢ 
-¢ 
~~ ¢ 
-¢ 
~¢v 
-w 
~~ 
a & 
x» & 


tee eF & & oe oe 


1 ¢@4¢44¢4644 8&8 8 wo wowowvwvwesedttitt@tete es 
V@eeaeae ee aewmreevewrewewveredttesb bse 
rtd ae 4 HR meme erewewvrwevveredteerteibaaa 


Pare Pew @weeae oe fb fF 


 ~ a awaaaeaeeaa4 6 4 


? 
q 





(27 = L1® C€E'0 = 9M Pure OF = ay ‘OG'Z = 414) 


0S'°0 = °U/Z 18 aDUR\SIG ASIAaASURI L SnsiaA XVN A /A €# YAS “ZT ainsry 


Cu/A 
Ov O°? 0°0 OT- OF- 





43 


(ZZ = Le CEE = 9M pur Of = dy ‘0G'Z = 14) 


09°0 = °y/Z je adue\sIGq asIdAsuel L SnsiaA XVYNA/M €# 49S “ST IINBLT 


u/A 
Ov 0'% 00 O2- O%- 





(ZZ = L3& €EC'0 = 9M PUP Of = Ay ‘0G'Z = 14) 


00'0 = Cy /A ye aouRysiq [edVIVA sNsIaA XVNA /A CH JS “6 BINBIY 


XVNA /A 
on 8) 0'0 GO0- 





45 


(ZZ = .L3€ CEE'0 = 9M Pur Op = ay ‘OS'Z = Id) 


00'0 = °y JK JE QOUPISIC [PONIIA SNSIDA XVNA /M CH YIS “07 ainsi] 


XV /M 
LO 0°0 L°O- 





46 


(ZZ = L1€ €€€E0 


aM PUe 0Z 


OY ‘OG'Z = AW) sautpweads FH aS 


Ou/A 


IZ ansry 





47 


(22 = Ld €€€'0 = 8M pure 0% = ay ‘0G'Z = 14) SIOIDaA AWDOIDA FH 1S ‘ZZ amnsiy 


aa @#@e tHe agaewseateees#sae@ 


aa naam aeaeaeaieteteeaeaeassteee 


aaneaanaaswsaeaeewey 


@ @4a 24842 @B @2 we weve 3s 


eo eo 4 Aa KAA KRM MH OD eT se sdsss 


i a i A a a a, a ee ee ke er | 


> PP hk ehh MRH OY eee ost Sb 


i a a a ee, a a 2 


| i A a A a a er er Ye 


| A A a a a a a a, ee ee Te ke Yr ae Ye 


yrPrPr rem PrP cP oH He Meo wdeysseabrkass 


yr PrP PrPPrrPrP?r ® 
> ry Pat ?e PF 
yr rr HTH 44eeee se 


-=vywy> db bh b b& 8 


vy > dss bb th 8 Sb 


>» ' abd db db 8 


&* a a & & &4 & 4 4 


> r, Fr wm & wm & A A 4S 4 


ovr FP Mmemanseea @ 4 4 


saab ' 8 6 
oss shb bb b&b & 


= aa 8 SS 
ovearer’ 
a leintaenl all ail ait aa A ee 


- 


x 
. \\N\Ne 7 
» CAN. 


A 


» 4 
a 4 
a # 
7 
f 

| 
i 
¢ § 
’ 4 


aaangnuanaa 
basen tuanaa 


(eA mana nmanm nw ea te 42 4A HG aoyn gb bh_ hb bbb 


a vrvreEewtenrewweeebhtbt4a4 


aetgcas32 @ wea SB a £ 4&4 4 @ 


aaseansunuauaa 
oaaésééssb&8 sb b 


ee a ee ee ee ee | 


avPrPPreF HS BHABHA HD 
Pe ee Pk fe ee 2 ee 


qveaeeerd ff 
Tee eeaee 0 


4é646464 6 & & 8&8 8 
¢e¢@#4é@4t4 46846 6 bb S 


@qa@¢€44444 & ROP we weeedet#e bb & bb bb 


q@*<18¢424844 4 & DP ewewee¢ ¢ ¢ £(4i4tl4lAtlk lt 


qt @¢€044 4&4 4 BD BRP vwEPwrPwrewtwttedte*ebtebBba 


£444 4&%4 & & HB wMOwPHPweweweweteltetlkh4lhlUC hel hkllf 


2404648644 & RS RB wmPeweewrvwveevetvr_t@+e_a@tes 


qt 4444 & S&S RBemerreewewevsveww@et0t1ietaa 


ye 444 & &R BerwrerevrPwrvsvse teeta 


a4 4402 & © BeBe vwrprvrwreewevwrrwtceveewiaéash 


a2 4 ¢@ @® ®@B Be we FP Pe wPereerewevrewrevd¢ 


&* axa @& ®&@ we earrnrererwreeeverrvreFe#Fe#eeg¢g 


anmhe Bn @ F BFP B&B FP PP PF PF F PP CF 





48 


(ZZ = L 18 ELE'0 = 9M Pure 0Z = AY ‘0G°Z = AH) 


0S°0 = CU/Z Je adUR]sSIG aSIAaASURIT SNSIAA XVN A /A HHIBS “ET aIn3siy 


Cu/A 
Or O'S? 0°O OG O'v— 





49 


(ZZ = L 18 €€€'0 = 9M pure 0Z = AY “0G'Z = 14) 


0S°0 = SU /Z ie BDURSIC| aSIBASUPIT SNSIBA XYNA /M FH IIS ‘PZ aINBIY 


°u/A 
Ov 0'% 00 O7- OF- 





50 


(ZZ = LI €€€E'0 = 9M Pur OZ = Ay “OG'Z = 34) 


C0iG= JK JO DDULYSIC] JBOTJIDA SNSIOA XVYN A /A FHYIS GZ aIn3iy 


XV /A 
G0 00 G’0- 





(ZZ = L 1? €€€E'0 = 8M Pur 0Z = AY ‘OG'Z = A) 
00'0 = OU/A 3k aouR\sIC] [BIIVIDA SNSIDA XYNA /M HH YIG “97 BINBIY 


XVINA/M 
LO 0°0 Ur 





(LZ = Li? €€E'0 = 89M Pure Oh = AY ‘O8'EL = 4a) SoulTWeans GHIeS ZZ aNBIY 


Cu/A 





o2 


(ZZ = LW €€E0 = 9M PUR OF = AY ‘OS'EL = 1) SIOPA AWDOTAA GHIAS “gz aIN31y 


SU/A 
O'v 0 0° 0'o— ° an 





°U/z 


(ZZ = L38 €€€E'0 = 9M PUP OF = OY ‘O8’EL = Jy) 


OG'Q = CY/Z Je VdULYSIC] AS.DASULIT SNSIAA XVINA /A GEIS °G6Z DANS 


Ou/A 
Ov 0% 00 O¢- OF- 





oe) 


(ZZ = L138 €€€'0 = 9M PUR OF = AY ‘O8ET = AH) 


0S'0 = OU /Z Je JdUL}SIG] aSJaASUBIT, sNSIAA XYNA /M G#ISS = “OE ainsiy 


Ou/A 
0" 0% Oy  OC— Cine 





56 


(ZZ = L 1 €€€E'0 = 29M Pure OF = Oy ‘OBE = 4H) 


00'0 = CU/A Je aduR\sIG JBONJaA SNSIBA XYNA /A GHIIS “LE BINSIY 


XVI A /A 
9°0 0°0 9°O0- 





a7 


(ZZ = .L 18 €€E'0 = 9M Pue OF = OY ‘O8'EL = JH) 


000 = OU /A Je aoUR]SIC JBOTIVA SNSIBA XVNA /M GHYIS ‘ZE BINS 


XVINA NA 
LO 0°0 L .o- 





58 


(27 = L3e €€E°0 = 9M PUR OZ = AY ‘OS'EL = Ay) SauTTUIRANS 94 IES =“ EE ANDI 





59 


(ZZ = L1® €€E0 = 9M PUe 0Z = AY ‘OB'EL = Ay) ssoyaA ADOOPaA 9419S “PE aN3Iy 


OU/A 
O'b O's 0'O Oo O'r— 


aS 
td 
a & 


4a 


_ es —e “_ “~\™ BW & 


\ 
\ 
\ 
\ 
| 
| 
b 


Ff” a sll 


ad 
9 





fV/z 


60 


(47 = L3® E€E'0 = 9M PUP 0Z = AY ‘08'EL = JH) 


0S'0 = CU/Z Ie DDUR}SIC| ABSIaASUBIT, SNSIDA XVN A /A 9H YIS “GE aINBLy 


Cu/A 
O'v 0% 00 O7- O'F- 





61 


(ZZ = L1& €€E'0 = 9M Pure 0Z = AY ‘OB'EL = A) 


‘VO 
0S'0 U/Z Je DdULISIC| aSIaASURIT, SNSIBA XVNA /M OF J ‘9¢ a1n3 
aS INO] 





62 


(ZZ = L3& CEE = 9M Pure OZ = OY ‘OSE = AH) 


00'0 = OW /A Je DDUL}SIC] [POTIVA SNSIOA XVN A /AQ# 48S “ZE ain3iy 


XVNA /A 





(4% = L318 C£€'0 = 9M Pur 0Z = 8A ‘08'EL = A) 
000 = Ou /A Je aoueysiq Jessa, snszaA XVWA /M O# IIS “BE an shy 


XVM 
G20 0'O GZ7'0- 





O'V 





(€CE0'0 = 89M pure CF 


= 2 0G'Z 


= 4) OL 


SU/A 
0'°O 


Lie saulweans 6¢ am3ry 


(€€0'0 = 9M Pur OF = AY ‘OG'Z = J) GL = Lie Soutpuvans ‘OF ain3ny 





(€c0'0 =a = ay ‘0S’ 
M pue Of = aw 0S°'Z = 14) OT = LI SaUTTWRANS “[pamnsiy 





67 





(C€0'0 = 9M Pur Of = ay ‘0G°Z = J) GZ = Lie Souyweans ‘zp amsiy 


OU/A 


4 


°y/zZ 


(€€0'0 = 9M Pur Of = ay ‘OG'Z = J) 97 = _L3e SauTfWLaNsS ‘Ep am3ry 


CUu/A 





OS 


°U/Z 


69 


(€£0'0 = 9M Pur OF = OY ‘OG'Z = A) ZZ = Lie souyureans pp ain3r1y 





70 


(C€0'0 = 9M PUe OF = ay 0G°Z = JH) 8Z = Lie SauTUVaNs ‘cp amy 


OuU/A 
0'O 





71 


(€60°0 = 9M PUR OF = AY ‘0G°Z = Wd) 67 = .LIe saulpwreays “gp amn3ry 


°U/A 
O'O 


y)) 


. 
v NSS 

es aor 

ano 


' e . 
5? @ ay’ “al . 
het ¢ fy Sen 
’ “ Si , 
* 5 E . 
ac 
« 





Hos 


O'V 





(€€0'0 


3M Pur 0F 


Ft ‘OGL 


= 41) O€ 


Le sautpueans ‘zp amsry 


O'Z— 0 


Ue 


REFERENCES 


Brackbill, J. U., Kothe, D. B., and Zemach, C., 1992, "A Continuum Method for 
Modeling Surface Tension," Journal of Computational Physics, Vol 100, pp. 
335-354. 


Kothe, D. B., Mjolsness, R. C., and Torrey, M. D., 1991, "RIPPLE: A Computer 
Program for Incompressible Flows with Free Surface," Los Alamos National 
Lab., LA-12007-MS, Los Alamos, NM. 


Neubert, D. E., Jr, 1992, Trailing Vortex/Free-Surface Interaction, MS Degree 
Thesis, Naval Postgraduate School, Monterey, CA. 


Ohring, S., and Lugt, H. J., 1991, “Interaction of a Viscous Vortex Pair with a 
Free Surface, " Journal of Fluid Mechanics, Vol. 227, pp. 47-70. 


Sarpkaya, T. and Henderson, D. O., Jr., 1984, Surface Disturbances Due to 
Trailing Vortices, Technical Report No. NPS-69-84-004, Naval Postgraduate 
School, Monterey, CA. 


Sarpkaya, T., and Suthon, P. B. R., 1991, “Interaction of a vortex couple with a 
free surface," Experiments in Fluids, Vol. 11, pp. 205-207. 


Tyvand, P. A., 1991, "Motion of a Vortex Near a Free Surface," Journal of 
Fluid Mechanics, Vol. 225, pp. 673-686. 


74 


INITIAL DISTRIBUTION LIST 


Defense Technical Information Center 
Cameron Station 
Alexandria, VA 22304-6145 


Librarian, Code 52 

Naval Postgraduate School 
411 Dyer Rd., Rm. 104 
Monterey, CA 93943-5101 


Department Chairman 

Mechanical Engineering Department, Code ME 
Naval Postgraduate School 

699 Dyer Rd, Rm. M3 

Monterey, CA 93943-5108 


Professor T. Sarpkaya 

Mechanical Engineering Department, Code ME-SL 
Naval Postgraduate School 

699 Dyer Rd, Rm. M2 

Monterey, CA 93943-5108 


Naval Engineering Curricular Office 
Code 34 

Naval Postgraduate School 

699 Dyer Rd, Rm. 220 

Monterey, CA 93943-5109 


Dr. Doug Kothe 

Fluid Dynamics Group T-3, MS-B216 
Los Alamos National Lab 

Los Alamos, NM 87545 


LT Craig F. Merrill, USN 
c/o J. F. Matthews 

109 Laurel Ridge Rd. 
Canton, NC 28716 


Vhs. 














Se 
OO 





DUDLEY Ns 
NAVAL P 


MONTEREy CA 95945. “3101 





“HOOL 


















Meta tstd kid re ihe i de ote oe 
fo bm are a z A ‘. 5 ’ ‘ Ste foe. 


bly + t-te Dit. ‘stebe Z : -. oan ty = 
ISLES ITU ated re RRR ee CM LTR DUDLEY KNOX LIBRARY 


ren Bets . b AP 4 f le 
A wet the Pe pasdg AS ae , he ‘+ ; Pa tate Ul. eA ster og: i Pike pe 2 Ate + - 1 ap ? ; ai 
pk S A aco Wem ae 4 i baie Loi Rae . La SR Cue CHOTA T TUT AUT TE TOP TS A 
i PGi: Asn gh F erral er - Piste Sy fies, ety x ire? aD eine iG ireenp siti: ron nt & a ae Oyo é ait | { || 
tab hee, belt t p Pali le be. GUNS Ort al eee hse 7 / it ait 
Le a ts a eRe ts tah.» x: o a +." - 73, bv atest, a bat tech t sof. ge eet % ‘ } } | } } | } 
ebete he ub eritat td & 4542p be OF 4374, 82 MIT: { ere pie oe hie Pcait Siceer ees Pre ity rede ee ay fo : HY | } | } | 
4 ¢ : oe 4 . . , Pa } | } 
ateb gett a tehce ses ky ine re i ee aybasenes is is : prt at we pen me Pod ot ait eee (OH 
eh, te 


SoM Ash hed oetaheasd asta 

eae ot rh 

fp daub A tte hE, A! 
nike Besos 

























" DY ae * en W.# « > . . 
Pa Ra OP AES ak LP pb behes © ar * ‘ « ar : * . P 
Rare se dus ht Oae 4 4§ Ps . * 
Crt: 439,00, 829, 9 ig? ae a . ae a 
4.4 cy a . b ni et had . ae) . . 
i he ya 8k rons Gnd ot 4 to a 5 : 
“ules Whee . * iT a 



























1 


oe 





Fhay bg : 
seapenveraidietsin rT r sit 
a 


gir iaaetcd 
at 






















































































cote . od 
4 er ee tr) & 
PG DBs ey Me tein we . 
irs astakersee rae) 8 sateiatth tee i : 
wot Bewibs 3 ta agit rf 4 ‘ 
hel ko rf i aaa . = 
Pm al s. 
hs hls} oe 4 a wae aatt e ae 
le %. 
ar Bad nog pe bisa pecie be i tale 3 ee eters 
tera Pv AS At bah 0. eb hia Li yaet 
cect thw, einen tale ta Cho ' aE *3 . ° ° 
ae Oak. a eter sae J ' 7 
Ve AL Fe ph ee 
tal pih-le eee asehi Bat r ane 
ati ry coud tg BhO, nity wie mckadh, fat. 
fare cake Mitte SOTO | ‘ resi «fp uah95 76 Anbes i . oe 
robin’ Sathengiok wales iy elke Wt ncbl hI ye ih aie vehi ee es is et Aces 
ce aad foe 8%! he fers yacht iw iui act ayh t. oie xi abs OT “ a * * 
ti! Tr a bh ebrtahe yk Ct bch hs 9908 & Pai ery st a Pa rs a sb ihe e. « + 
(oe WF, iit Es Ft Rat. ww Ty, ig ag t Bg hele et 9 ea dstaratole “tight ta? ben gt ae e 
Ase tanp hey: ih LcAyn es Palghs debs idm itele Secor aie i abe y shout Pe wr Te Haar thle} HF tenia 4 bit - 2 
eka Sh a Pipe aay jist each Bie bem “hy by iat boedh beens: are! rigs fi whee; abt © , es 
debs Vat na4.s 4gneh te spakiia Hite viayiisi hes add sce Ps 
Gyre rhs ES i, Pay ae eat 19 Aven eile ie Wenge CAMs igs fig egtveg’ we - 
Pree ae oc eer #4 tebecicet ty ee fps # z elaqiiesine eae, 





Pia wrth ie 


teks red ee Heed Abbey *, ah ates igh hee Seb a 


* A AL ck +s 
AT ACF y ‘ fl el K y, eahenaled 
aR ait sities eine ern F ha rapes 5 ue nya it hie ‘fice! bape seiphe Xi * eb ef wha . Sees 
ri 


% 






























hata hal sdoiatatsnss chats ea ture: Py il * Bae ahs al i; wishes MAL 1s 3 e 
i . oi tie “iyi g che ry 
A path contin Ne aS ae Noi sents a a eat Ss ene 
had ey 4.9.8 hss ye tate ts R. abe ree, ans ry Ret oh me cer hod. ced she 
pine te Fionn me ii fa te b> Meer 4 pte . 
a} ier 4S Chel - girlaiel 





tot ‘i 
Rar hah inh ieee i is 
ry ey etek ret Pre be oper Lt br} 








sa} arene 
iy he Phin aisha or 
A) aah rahe Ay Pe ae he Hetty. af; 




























































































































































tte b 
i : * ; 
2 is bet "a at? be ia in 69 2 a : 
fea ins inig? aot aoe 4 4 Sita ge aN vita yer ryt asebee tae : 
noice Nee are Oe 48 rT be pe : ih Me citer rath tig’ i c 3 
he i Tae am is 1a - sa a . 
VE ee Ae A within ke ect fe Lar Re sean ATR wad fe on nS 
aie tate ine iets esha * preci ob dae afi ? ss a) ac 3%, nh - ‘ 
be a ts heh seth rake lat Fa iey ie wy t a rs * 
titans °K Be J ciples, aers &, ( au “é Mie ha AES ids o 
ar i ietebs ee wane Sd hie yay a hasty Laer 
: ch ae a RE 1 ‘ 8 san! isan i pa 
ehse. ¥ let ae 7 Pd . 
tes te alex: be 2 ae ang batts cial ne setae Pe ys . 4 
* on . a rary FE ed Ht *.8 
eva! grat iat aa rite ae es nate at Me ies Rid ai ta : 
4 Tit et ea, tg tho uptake Ayia tats Le i re yi 435 35 ou) at iv } ts 
af dae fea ge Mg ST aie eis PF rdghat #30 By, aby ; - $ Vet ate ws TY so) Ye oe 8 
t iiltg ‘ iy ‘, SAGAS AIDS 11 5® *§s yee Soe a 14, ? 3a" 5 “hoes ot. *3 “He ren 2 «tar gee 
nt ape ee aed ea LEite Pf es OCS Be TY Se a Aig Gi ‘ ST a 
Bs hy UM Raat Hite Sis at Hess pi te Per ee - P an i = 
trtyb tis. aay By A 4 Pihely ek <2 pe ah ey Sytagoty ‘5 a egies ee ‘: fe) DLP ooh y i” as AOR ee 






gneh 0.6 GOS Or ee 

Re tis nee ue wl bape § os ol m4 

a wei vi Tlyartgapedss a t.F ton 
PEAT soe Fhe f aarct aaa 











o Be 18? ace thtee F pt hent “a ye ¢ a ne 2 
sy tt stapetel @ fas ¢ etack- as 2 +3 anf 





: fig a a) 
ae él 


a vn 


La ae, om © 
pase , 


x dens aa EE 
is 


BA hs bake he 




















, 
Ca ae PP a tienen ae hee = 
. 






















































































ei ry tte RARE Ape «bse Ue » OA | ie 
aptaiat Lae} ay tbe 

2 uae Bhseits? ie % i co ‘ A a 
eee Ann ae eH csteig) F dete ‘Siz eh th init A Febdg ior tab ahi weit’ = 
Oe 7, Dabols Wy bg igh ae seh ee pik eth Po baa 35 
ere etd Lid whi seetes ion aoa tae Hoth rhe ftiikid mga t om rae ky by ni ye 
PP yg Fins seine Peres pl hate My hee DoF Hi eri ats veh a hs E saet qidect b i: ’ 
praia ae aie Sahih, oehen. ety 1%! gesd hee el Le 3 “aye e 5 
hai bene af Ta eer ph “ >) eat i! Ge Pare Sry ar) Pyeyerert seth i 
HAs Testo sana % rads H th jie i pies * pee ee os fy 
Sema ¢ 2 s pe tet hg? eit ; Ai iin 
ghey SaPod Mh ioe Witgd bet) yi 4 Behaeakat rte 


lobe org PS ie vias 

pabeineanse nyse 
Fp LM oy Aas 

ip 8% ate 

Fi id ue pest 

eee 1S 
Ea, yates fee 


ee taf Rate 








i fs 

SFE34 ite ite % ‘ rp! ? is ‘ “ Sats Feith : ek 
RSA ies ae ER ee 
ae % oa ’ hh) ws * i 3! ‘ oY rit iHst iit 





































a 






































































eee hte uit bo ee ; 
i e. a . ee 

retet deh ei teh! ri 2} = es fateeta age ek ple # * 4 seeped agers aun e «° - 

fs 4 5 q Cees ea « * oe hens oe “ a . a 

Abin ta oes, i Pell is af Mi ey sh ge . ct a bpt + emr g"Be & ste ie os ie at ry Pi ee . 

AY ay At een re eT oe ir eer a | a 46a ee ea) oe 6 om a) ee 
phe jit iit ¢ 20 of «pew ato fC Py or. we oe = = 4 se 

sient + wth arate ae . ve . 2e* oy « wae eee a . 

‘4 venting? 8h th0.8 « eo 86 a "8 406 fea ¥hs ory - . - « as «ae * * ° 

als if sah ahs iy hea ieeghilae Pint acme 4 te Dee 26 eth = £ Bee Pere nk ave 8 ea @aee Si ee a é 

Py aes aay : ; b: pribete LA ATA Soh s Bead ve aie = Sistes@e  lagele 8's aNs ‘S : s = 
dias Pte aEnPal ieee ‘ HD fescipritys | Nici } ans e see ce aa Naietcaare. vo Fs = ors 
wi atie Fpaeak: Pads ata ds 545 pi ghee oP tabs rs ‘ 6 ©. ee he oe . . 7% 
Tang! ites [rade d ee Vibe rk Hee teataws tt P palegy nto ee * s Chr ets ee ° P 
4 ti eeeiten Fi tye, ih ifs! et Teh oO Gee ees Pyar aaatite oes we tate 
oe els wares j o Bek wan Uy oa a Sip e 
dca rye One Pars ayrary bd > "ate o 
cong says she! Pa, eNeidtrah 0 9.4 ip Pe eT ee 
ie agg ya rh aoe ht +E woe, { ee. ig 2,04 a * ue ' 
pa i Pyige pateescsss fal ‘ actre J 
~ he 2 , be ar 74 4 
5 ate eed wit; 4 : od oka 8 « 
atts Syed e tes 
‘ 










eit 







Ui Hite Pi 





daisy , ae ie 







=< 


















te 
se 
hey 
ae sf auaee 
yu ”. mY Bageze af! 
3! 
5 Pes ; ey a 


fie cat oi 


toe, 











é 
: wine Hie stants 
Fe ppv. stot if ite 
Supe Sa MT Be 


ales BL Bente 
ae 








B Et 
a 







RERDGE4 PNG Potet a> 
pevre sets tr et Ha 
eee? BEA i is 















eae 
wis 







































te I 
e 35% Patt, he 
* ita GF nm Ue ae ee t 
e557 sod aha: ass iedys ’ neeichs | Bs <4. 
3h s Pyayhey Ste aed Sitestt 4 
ary iH Vv? i? i 
Vy sary Fyiyh aby Pyema tua 
Grae oye Bed sey Pytar Fae 








J wee 
Sed vy 
& 


ay ven} ane ae o cae 









aes b, 





















































‘ AP % ye £ 
nay ili sa mob, cae santas “we aa z apie, ¢ x ‘ “ 4 tf 
; man sates f a : i ae Kea uns b tind ; 
vst) ne i PT RE NEHA Men Mete HEM AEE dpe thi Oat 
A Renan y eis neat chy tage yim! fryers Teleten ys REL bs ehies de wae 
Dena? a aa an ce ee 
r a ee rk ' 











an : un 





ee 
fa bight oy 7 ae 1 












pet 
tal 
ai : 






















rey} 
es 













































at ii 





te! yey" Ht 5 
i rab eres r Cea i) pe a i reer 
ital “th oe 7 REHE Le a Ht en or Hay RAN NSE cad OS 
= a Se Ae Tie! " i acfaaes bpyf om iit gryrets Wah hs: Pith at Oa erate red 
ere ee ae a yt oh ihaee ore os iy y re Mytde gt i ; aed 
in a \\ yey vii it y pele ny iy fk A rice pielyee os 1 Ly gt yes fdyre ia 4.3 hese 
2 : hy hal sty 4 5 shy tg ‘ayy eras Ped gh i aes $ ait Hot iat emer amety f) 
i Al ar ann ne cae ah iif fs Setsdys ay ree! gh yf na say ae ays Hy ue: wi set 
5 : Tat i) ‘iat 


Tete (ry 


ES res 





pases 







cat tctnlsts an 















































vay ty Cie 


Matt 





























it ty aa i id ty te i! uM seats BURA : 
fy et i fe n oe cai ore sf ie oe eee ait Avy Aon Rie ey en a ey % ie ey ro 
ihe a Peta pent Ne HH ? yyy ae : sy yew a 

nets : hepa | By y 3 tn nit ifs) ebely sf ht aii oH ae Seu styes a is wie er ‘ bite tog his a s8 
J i dabei cfate ae mec me ithe AX meat Bites thy yt fetta yh ay dena A) MBIA 4 ‘ it st! Neat « sradge ; 

aie xe pe ye iy a mF 2 ata P58 4 Sie! } MD Nis ein thd Hat tad mitt ART: ry 5 Bet Hy 
a mie uh 3, Seah ny oe Rae nl w iyi “8 ie? parity v4 ane welt sty ars 
vytsky yor ra stent, Ae: Ree ‘en ni eri iste sh ve ¥ 


ay nie Hit Nie Kin fais! a te Mh Dia eh ith nin Hes te 









«4 12 ben eten be ‘ie sea if fe ry He ¥y A sat 5 aj att. By eNghe € 
sens Ae Haat i ' Me wy ny plat tT uy Hite) ae is dl ey oy a iy set Re it hh nay 
yhhryn i fan nic yur s¥5 ‘ 







iy be 
a nhesh rye J 
Ht We Pay: 


y , ote *&: 
: iste =A Dee au Haye : 
Ses La whe it 







on RARE: Ly} aha EA, UM tad | F Merey 13 itd 
Boa ne aa yates te i RA 


Re Tan ease Tate 
peek! ae Heh Basia hale 






- rh tt 
aie 









a . tye Ui isi: Woy’ J ats : Mens if * 
Sat any ie ribet fi ay tie pune aes teste ays sana See fat a ci fondly, Apne) eae shat, 
PMA RM AHA mn tela et see aN CU tetesytaty ote perc teat are, Oe 
i hl ia Leet th Dat ‘nA pyataisenert flee ayn ie Lats Lao c ras ns ese Han Rana ook ames eatin tk 
; : pi! eee y te nasialatanet te iu: eg leat edt sy Awe comes UEC, 


tiftibes 


firttlnygae AAT sia UTR Peet Nte *! Vitel ah ahe th ws z ai st pk a a 


thew, UF 
ran 





5 a 
MA toy i Mii ia “iy 


















tidus ry oforah mk wie Gta pe © ‘ eeeao kya © 
ee) ifety Se i noe f ytd Ene ae 4 anh ne My: RY bn tat BS Hels wey ttt Ak 5 Vit: hays sb "4 tay fy ; my sy 5 jt! 54 J tyhe Lit . an 
if ty i} as se ! rh Veo" wt t ghee rs alysetet yt ve es me a ital . sh ae ADS Paseet Bee th vee ii Le a ; =. rr *! “aL ' fetes oh ets . 
Nps whan? yeaa nit *\: “ie a A a Wayige Asay sate DUH ty THAR rns yfeferds ytyin's shhgretgten ts uly gtel fed gil’ y et 
Me Rt ih ia xh i ity gipatel ket “st ATE sie hela f \ Fo 3 in yaa y - Ae te iM rseaiy yest! al Neh, n hy ten stat oy * 4 fe)! rbees utte"s a Potccs ey « vv seas con Ly + “he ie 
z i Yi sit a febey iy? Js yt Hts ig Ea oo ie fe ae rity Wty Y of oh, vod tara]. Kiyed Da pha pes® (hy os s?y ad © ae wna @e ry aus 
htt mi Pa aie a a Gis i MG 1,7 ty vans Hees nt ‘<* agile Ve he 1s yet mis reds a%po , q dt« + Ot ero At | tt Bey ae 
r nan ty camer Map ie ; ase ayes tie ray yay sets ah Mets ety tat \ vty: Jyld alate? 7.8 hie ae tie ARE Tee ; 
an gees righ sit by dike 14 hind Wy re fe nat v Ha, fied fi" shyly a} oe 4; nt yes wy a ay! vials ’ ‘ie NSiterse re sek earl 4 ue gts by gia sSs. $08 toys x Re Sie | “a! ve 
pen thrits tol bie yt mi Me i Ve ect th ey Hi! me at; sont aR abetabil Weta mae erat hed He Hille Nee tee ONS Madea RY Sk BE., SO ye 
Jee eed ten) mea yt te) aa ite iy Bh ir ean Waite “fy YY abyetetan te { “the rie we siytaie ‘ afte viol pighus se = he oN et felt Fe A : iy ae 4 4 . ae eb ‘ fae 7" 
iy" shy vi) ate : iy ? et mani inane re, it eR aay tte ere! ns! ‘ied i ‘t Abe? «8s! tty us ase chen at gtnig ede o & Wee std at mat ree 4 ray mt x33," i a ee 
alata Be Pert yi pa yatah ie Pete yee Niet ial rary ! M ait evant: k roy Mt; er ‘ “erent htt eiud ov cayhbhege splerel 6 ORCA ERTL 0 ds A” oper mS 0 4 ce 
Hany nytt ale iy wy cit ioe ie i hay hae te His a3 Si Aaa | ‘ ate pee is ear LEN REG LA ae | k Mery aoe ne AL ——T- . $ 
esi Hs rt! i Heit hyty : a Hate CaN ie a ey ete heen ars afyisut my dean tiniie He “: abst thks . 
cash ‘ ra ¢ if HP feat PEERS te oie ae Kies tp) at Nat, pens a eres 
elytgemnths wry Hecate iit Mh it HM i Hrd toe Pt It eae OOK ae ar abt 
mn 4) UH sheila iat hy Vo. reg faded tug bag sk ae : 
ie wore in Bs ny ; she eat s A Nii fein s Byrne! wy per abe sky ms Ji ath v: 
yh ait 


sah 





Bist é PASH ai sis i Pacey A tc) thts fing et pests 
say Ay Cit my ate “ni U HATS belie lag a 3 thes sles #51) He | ROY fh iy a ihe Pa 
ee iat ast rae Ki iy hi sae fon Pe teat : 


as nits tit, HY vaatytt el 
Hts mile ye} pa by! i hing ef | 
- Mths ny tity nt \ Hh i a yu the te 1°s prises oat 
ee fa Bis ee cnt ante ve RRR OCR 
shy ses nn wi ile Nees ries 












, ae 





ely ty ai 


Pray 
whys ie ey ae 5 hat 
by: ai rs aay i ve! 
ry 






i“ hey 1 ‘ts it t 7: 
,, My 8 yess te tal Ries : , Wh: " in ive fee ae u HH 7 Pa Ye) 
ie ROME ee bites eth ; it ; Wye 
és Pe tie aS wa ad GR sh ine ba cae Noes 
of oF eps 


Y ‘ MA) A) sr i it ay Pree - ‘tot of rit A qe ot: By seypha re > fy fi x 

cbse Dae ahet te A wie siya a PE RR Moa 

ai dais ie ee rier ee gs ame 

p vey ha 1 HA ‘ ane r te ne! oy iat ie fps h its i f sry tatal 
nue sha i i fu iy we Atay aus ua it me vit Hite eh \ t bietiate 


chy megane Hide, FUN nue ptt NE "3 


Sate Mis ph § pan av . iy H tH, 
Cy ink ye rat ytah 19 fae vat v1} i ve ity ty} Pe tes a We weary pene e furl if , 
fe! arise ‘ steely A Tye oF ae sn Paes bias sh ny") pi YE hes yea ACEH wi eh fey ele! 





